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'_- ■ The goal of the present paper is twofold. The first target is to present an enhanced group 

classification for variable coefficient generalized Burgers equations u t + a(u n ) x — g{t)u xx . In 
other words we complete the results existing in the literature. The second target is to show how 
simple is the procedure to solve boundary-value problems using Lie symmetries. 
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1 Introduction 

The Burgers equation, u t + uu x + u xx = 0, is one of the simplest nonlinear (1 + 1) evolution equa- 
tions. Nevertheless it has a long history as it was already known to Forsyth [11] and discussed 
by Bateman not many years later [2]. However, it was a serious contribution made by Burgers 
which led to its present name [8]. Burgers equation has been used to describe many processes 
| in fluid mechanics and a variety of other fields which do seem to be rather disparate. It can be 

■ transformed to the standard heat equation by means of the Hopf-Cole transformation [9lll3j. 

In 1989 Hammerton and Crighton [12] derived the generalized Burgers equation describing 
the propagation of weakly nonlinear acoustic waves under the influence of geometrical spreading 
and thermoviscous diffusion, which in non-dimensional variables can be reduced to the form 

ut + uu x = g{t)u xx . 



In this paper we investigate Lie symmetries of a generalized Burgers equation of the form 

u t + a(u n ) x = g{t)u xx , (1) 

where a is a nonzero constant, g is an arbitrary smooth nonvanishing function and n ^ 0, 1. The 
enhanced group classification for this class is achieved and comes to complete the results that 
exist in the literature |10[I30|. This generalization of Burgers equation was recently treated by 
Abd-el-Malek and Helal [T] who used what can only be described as a very ingenious method to 
obtain the solution. Here we present an alternate approach to solving the equation with associ- 
ated boundary conditions and we are of the opinion that this approach is easier to implement 
and somewhat more transparent. 

The paper is organized as follows. In Section 2 we determine the equivalence group. In 
Section 3 we calculate the Lie point symmetries of (pQ) and in Section 4 we demonstrate the 
solution of the boundary-value problem. 
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2 Equivalence transformations 

The first step of a group classification problem is to find the equivalence group admitted by 
a given class [22J. We have studied all admissible [27] (form-preserving [15] ) transformations 
that link equations from class (pQ). The results of the study are given in the following statements. 



Theorem 1. The usual equivalence group G~ of class ([I]) comprises the transformations 

$3 ri-n ~ ^3 2 



63 <^3 2 

t = 5\t + §2, x = 63X + 64, u = 6§u, a = —5\~ n a, g = -^—g, n = n, 



where 6j, j = 1, . . . , 5, are arbitrary constants with 616365 7^ 0. 

It appears that for n = 2 class (JJ admits a nontrivial conditional equivalence group which 
is wider than G~. 

Theorem 2. The generalized equivalence group G~ =2 of the class, 

ut + a{u 2 ) x = g(t)u xx , (2) 

comprises the transformations 

~ at + /3 _ KX + ^t + fiQ a 

t = ^ T j, x= it + 6 , u = 2q(q(5 _ ^ (2a«(7t + S)u - K1 x + Ml g - M0 7) , 

a = - and g = — —g, 

a ao — P7 

where a, f3,j,5, k, fix, Ho,cr are constants defined up to a nonzero multiplier, a6 — /3j 7^ and 

KG / 0. 

Proof. Suppose that an equation from class flTJ is connected with an equation 

u~ t + a{u n ) x = g(t)u xx (3) 

from the same class by a point transformation 

t = T(t,x,u), x = X(t,x,u), u = U(t,x,u), 

where \d{T,X,U)/d(t,x,u)\ 7^ 0. It is known that for evolution equations we have the restric- 
tions, T x = T u = 0, on the general form of admissible transformation [15] and moreover for 
the equations of the form u t = F(t,x,u)u xx + G(t, x,u,u x ) we necessarily have the condition 
X u = [26]. Therefore it is enough to consider a transformation, T, of the form 

t = T(t), x = X(t, x), u = U(t, x, u), 

where TtX x U u 7^ 0. After we change the variables in ([3]), we obtain an equation in the variables 
without tildes. It should be an identity on the manifold C determined by (JJ in the second-order 
jet space J 2 with the independent variables (t, x) and the dependent variable u. To involve the 
constraint between variables of J 2 on the manifold C, we substitute the expression of Ut implied 
by equation JTJ) . The splitting of this identity with respect to the derivatives u xx and u x implies 
the determining equations for the functions T, X and U 

U uu = 0, (4) 
gT t - gXl = 0, (5) 

X t U x -X x U t + gT t (^J - ahTtU^U, = and (6) 

n— 1 ~ ~ Tt TT fi~l 1 r,~ U xu X xx Xf 

anu - an xf +29 x!u:- 9Tt jq + x- x =0 - (7) 
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Equations (HJ) and © imply that 



U = rj 1 (t,x)u + T]°(t,x), X = ip(t)x + tp(t) and g = —g. 

Tt 

Here the functions rf(t,x), i = 1,2, <p(t), and ip(t) are arbitrary smooth functions of their 
arguments and r^Lp 7^ 0. 

When we use the differential consequences of the fourth equation with respect to u, we get 
that the arbitrary element n is invariant under the action of a point transformation, i.e., n = n. 
Also we obtain that rp = V n 7^ 2. After we substitute the expressions for U, X and g 
into the third and fourth determining equations, we can split them with respect to u. Further 
consideration varies depending upon whether n 7^ 2 or n = 2. 

I. If n 7^ 2, then splitting of (|6j) and (|7|) results in the equations 

vl = Q, v 1 + + 2<pgril = 0, 
Wt 1 = vl (<Ptx + ik) + wvlx and 

The general solution of this system is given by 
T = 5it + 5 2 , tp = S 3 , i> = 64, ^ 1 = ^5, 

^3 1 — ^3 ^ 

where 6j, j = 1, . . . , 5, are arbitrary constants with <5i ($3 ($5 7^ 0. Then a = 7-^5 "a and g = —^—g- 
The statement of Theorem 1 is proved. 

II. If n = 2, then splitting of equations ([6]) and ([7]) leads to the system 

ril = 0, r? 1 a; + ^) + lipgri - 2~aT t rfv l = 0, 
<PVt = vliVtX + tpt) + ¥>gvL - 2aT t (7]°r] l ) x , 
<PVt = Vxiftx + ipt) + fgVxx ~ 2aT t T]°r] x and 
(pa = ar^Tt. 

From this system we initially obtain forms of rj 1 and rf as 

1 <P a r 1 1 \ ~ a 

T t 2aT t a 



where a is a nonzero constant. The remaining equations for the functions T, ip and ip are 
Tt), " \T t 



*\ =0, =0, (t\ =0, 



the general solution of which can be written as 

act + P k mt + hq 

7* + 71 + 7* + 

where a, /3, 7, 5, k, //1 and fio are constants defined up to a nonzero multiplier, aS — f3j 7^ 
and When we substitute the functions T, ip and tp into the formulas for r] 1 and r/°, we 

get exactly the statement of Theorem 2. □ 

Note that the group GJT was found previously in [U] (see also [231I23] ) in course of the study 
of form-preserving (admissible) transformations for the class of generalized Burgers equations, 

u t + uu x + f{t, x)u xx = 0. 
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3 Lie symmetries 



We perform the group classification of class ([T|) within the framework of the classical Lie ap- 
proach |2ip22| . Naturally one performs the group classification for class (pQ) up to G~ -equivalence 
and for its Subclass ([2]) up to G*^ -equivalence. 

We search for operators of the form T = r(f, x, u)dt + £(£, x, u)d x + 6(t, x, u)d u which generate 
one-parameter groups of point-symmetry transformations of an equation from class (JI]). Any 
such vector field, V, satisfies the infinitesimal invariance criterion, i.e., the action of the second 
prolongation, T^ 2 \ of the operator T on equation ([I]) results in the conditions being an identity 
for all solutions of this equation. Namely, we require that 

T^{u t + anu n - l u x - g(t)u xx } = (8) 

identically, modulo equation ([1]). 

The criterion of infinitesimal invariance implies that 

r = r(i), Z = €(t,x), 9 = 9 1 {t,x)u + 9°(t,x) 1 

where r, £, 6 1 and 6° are arbitrary smooth functions of their variables. The remaining deter- 
mining equations have the form 

= (gr)t, (9) 
an9 l x u n+l + an9° x u n + {9} - g9\ x )u 2 + (9° - g9 xx )u = 0, (10) 
an( Tt - £, + (n - l)9 l )u n+1 + an(n - l)9°u n + (g£ xx - 2g9\ - i t )u 2 = 0. (11) 

It is easy to see from ([9]) that £ xx = 0. The second and the third equations can be split with 
respect to different powers of u. Special cases of splitting arise if n = 0, 1, 2. If n = or n = 1, 
equations ([T]) are linear and are excluded from consideration (Lie symmetries of second-order 
linear differential equations in two dimensions were studied over the century ago by S. Lie |18j.) 
Therefore we investigate two cases, n / 2 and n = 2, separately. 

I. If 7i 7^ 2, then we immediately get that 9 1 = cq, where cq is an arbitrary constant, 9° = 0. 

r = cit + c 2 , £ = (ci + (n - l)c )x + c 3 , 9 = cqu. 

The classifying equation on g has the form 

(c 1 t + c 2 )gt = (c 1 + 2(n-l)c )g. (12) 

Further consideration is performed using the method of furcate split suggested in |26| . For any 
operator T from maximal Lie invariance algebra J 4 max equation (|12p gives some equations on g 
of the general form 

(pt + q)g t = sg, (13) 

where p,q,s = const. In general for all operators from A ma,x the number k of such independent 
equations is no greater than 2 otherwise they form an incompatible system on g. There exist 
three inequivalent cases for the value of k given by k = 0, k = 1 and k = 2. 

If k = 0, then (|13j) is identically zero and c\ = c 2 = cq = 0. So, if g is arbitrary, we obtain 
that the kernel of maximal Lie invariance algebras of equations from (pQ) is the one-dimensional 
algebra (d x ). 

If k = 1, then g £ {ee l ,et p } mod G~, where e = ±1 and p / 0. In the exponential case, 
g = ee t , we have 

r = 2(n - l)c Q d t + ((n - l)c x + 03)^ + c ud u . 
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If g = et p and p 7^ 0, then 

T = dtd t + {\{p+ l)c\x + c 3 ) d x + \ ^ciud u . 

In both these cases the maximal Lie-invariance algebras are two-dimensional with basis operators 
presented in cases 2 and 3 of Table 1. 

If k = 2, g = 1 mod G~. The infinitesimal operator takes the form 

T = (c x t + c 2 )d t + [\c\x + c 3 ) d x - 2 (n-i) Cludu - 

Therefore we have proven that, if g is a constant, the maximal Lie invariance algebra of ([I]) with 
n 7^ 2 is three-dimensional spanned by operators presented in case 4 of Table 1. 
II. If n = 2, then splitting of (|10p and (jlip results in the system 

01 = 0, 2a0° + t 1 = O, t o -<7C = O, 
rt - Z x + e 1 = 0, 2a9° -6 = 0. 

The general solution of this system is 

r = c 2 £ 2 -I- cit + c , £, = (c.2t + ^ci + c 5 ^ x + 2ac 3 t + c 4 , 

1 1 

= -c 2 t - -ci + c 5 , 6»° = — C2X + C3, 

where Cj, i = 0, . . . , 5, are arbitrary constants. The classifying equation ([9]) takes the form 

(c 2 t 2 + c 1 t + c )g t = 2c 5 g. (14) 

For any operator T from the maximal Lie invariance algebra ^4 max equation (]14p gives some 
equations for g of the general form 

(pt 2 + qt + r)g t = sg, (15) 

where p, q,r,s = const. As in the previous case the number, k, of such independent equations 
is no greater than 2 otherwise they form an incompatible system on g(t). So three inequivalent 
cases for the value of k should be considered, namely k = 0, k = 1 and k = 2. 

If k = 0, then (fT4"|) is identically zero and Co = c\ = C2 = C5 = 0. So, if g(i) is arbitrary, 
we obtain that the kernel of the maximal Lie invariance algebras of equations from ([2]) is the 
two-dimensional algebra (d x , 2atd x + d u ). 

If k = 1, the following statement is true. 

Lemma 1. Up to G~ =2 - equivalence the parameter quadruple (p,q,r,s) can be assumed to belong 
to the set 

{(0,l,0,s'), (0,0,1,1), (0,0,1,5')}, 

where s' is a nonzero constant. 

Proof. Combined with multiplication by a nonzero constant, each transformation from the equiv- 
alence group Gn=2 can De extended to the coefficient quadruple of equation (fT5j) in the following 

way 

p = u(p5 2 — q^/5 + rj 2 ), q = u(—2p/35 + q(a5 + /3j) — 2ra r y), 
f = is(p(3 2 — qa/3 + ra 2 ), s = usA. 
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Here A = /?7 — a5 and v is an arbitrary nonzero constant. 

There are only three GjT-inequivalent values of the triple (p, q, r) depending upon the sign of 
D = q 2 — Apr, 

(0,1,0) if D>0, (0,0,1) if D = and (1,0,1) if D < 0. 

Indeed, if D > 0, then there exist two linearly independent pairs (5, 7) and (a, 0) such that 
p5 2 — qjS + rj 2 = and pf3 2 — qa(3 + ra 2 = 0. For these values of the constants, a, f3, 7 and S, 
we have p = r = 0. The coefficients q and s are necessarily nonzero so that one of them can be 
scaled to 1 using multiplication by an appropriate value of v. It is convenient to set q = 1. 

In the case D = 0we choose values of a, f3, 7 and 5 for which p5 2 — qjd + r^ 2 and the pair 
(5, 7) is not proportional to the pair (a, f3). Then we obtain that p = and q = v(3(qj — 2pS) + 
va(5q — 2r~/) = 0. Appropriate choice of the pair (a, j3) allows us to set r = 1. Then the residual 
constant s can be scaled to one by choice of v. 

If D < 0, we have pr / and can set p > 0. We always can make p = r = 1 and q = 0, e.g., 
this gauge can be set by the transformation 

2pt + q — y Apr — q 2 

t , • 

2pt + q + \J Apr — q 2 

In this case the constant s cannot also be scaled. □ 

Therefore, up to G^-equivalence, we have three cases of g which admit extension of Lie 
symmetry algebra by one basis operator. These are cases 6-8 of Table 1. 

If k = 2, and g = 1 mod G2 we get a five-dimensional Lie symmetry algebra which is 
sl(2, R) S A 2 . 

The arbitrary constant element a does not affect the results of the group classification problem 
and can be scaled to any fixed nonzero value by the transformations from G% ■ It is convenient 
to perform the gauge a = 1/n. 



Table 1. Group classification of the class u t + a(u n ) x — g(t)u xx , n^0,l. 



no. 


n 


9 


Basis operators of ^4 max 


1 


^2 


V 


d x 


2 


^ 2 


et p 


d x , 2td t + (p+ l)xd x + P ~ l ud u 

n — 1 


3 


^2 


ee* 


d x , 2d t + xd x + -^jud u 


4 


^2 


1 


d x , d t , 2td t + xd x - -^judu 


5 


2 


V 


d x , td x + d u 


6 


2 


etP 


d x , td x + d u , 2td t + (p+ l)xd x + (p - l)ud u 


7 


2 


ee* 


d x , td x + d u , 2d t + xd x + ud u 


8 


2 


g.g2/3arctant 


d x , td x + d u , if + l)dt + (t + p)xd x + (x + (p - t)u)d u 


9 


2 


1 


d x , td x + d u , d t , 2td t + xd x - ud u , t 2 d t + txd x + (x - tu)d u 



Here e = ±1 mod and p is a nonzero constant. In all cases a = 1/n mod G^. 
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Note that the group classification of class ([2]) was carried out independently in [10J and [30 . 
The results were presented there without gauging of the obtained forms of g by equivalence 
transformations and the case g(t) = exp ^2p arctan "^f 1 with 7^0 was missed. This case can 

be reduced to case 8 of Table 1 by an equivalence transformation from G%. 



4 Solution of a boundary- value problem using Lie symmetries 

There exist several approaches exploiting Lie symmetries in reduction of boundary- value prob- 
lems (BVPs) for PDEs to those for ODEs. The classical technique is to require that both equa- 
tion and boundary conditions are left invariant under the one-parameter Lie group of infinitesi- 
mal transformations. Of course the infinitesimal approach is usually applied (see, e.g., [3 Section 
4.4]). The first works in this direction appeared in the late sixties (see, e.g., 1, (>ll 19lf2"0"1). The 
authors of [HE] used the classical approach, namely at first symmetries of a PDE were derived 
and then boundary conditions were checked whether they are also invariant under the action of 
the generators of symmetry found. In the case of a positive answer the BVP for the PDE was 
reduced to a BVP for an ODE. Using this technique a number of boundary-value problems were 
solved (see, e.g., [H[28l[29] ) . 

In [3l[25] it was mentioned that for most applications it is not natural to require that appro- 
priate symmetries of a system of DEs preserve a particular BVP for this system but it is enough 
to impose that these symmetries map BVPs from a certain class of such problems to each other. 
Hence the induction of well-defined equivalence transformations on a properly chosen class of 
physically relevant BVPs can serve as a criterion for selecting symmetries to be taken into ac- 
count, e.g., in the course of invariant parameterization or discretization of the system under 
consideration. 

The method suggested in [20] uses specific one-parameter Lie groups of transformations of 
the independent and dependent variables of the PDE system as well as of all arbitrary elements 
which appear in the equations under study and in initial and boundary conditions. Namely, only 
the groups of scalings and translations are considered which can lead to self-similar or travelling- 
wave solutions only. After the admitted Lie group of scalings and/or translations is specified, 
the complete set of absolute invariants has to be found. Then a boundary- value problem for the 
PDE system is reduced to similar but simpler problem for the ODE system. Such an approach 
was applied to a number of engineering problems (see, e.g., [1] and references therein). There 
is also the approach in which the group classification of the PDE system and associated initial 
and boundary conditions is carried out simultaneously (see, e.g., |16pi7j and references therein). 

Initial- and boundary-value problems for certain classes of nonlinear PDEs were solved in pQ 
using the method proposed in |20j . In this Section we demonstrate that the classical approach 
is much easier, especially taking into account that problems of group classification are solved 
already for wide classes of nonlinear PDEs. 

We consider the initial- and boundary- value problem 

ut + a(u n ) x = g(t)u xx , — oo < x < +oo, t > 0, 

u(0, x) = 0, — oo < x < +oo, (16) 
u(t,0) =q(t), t>0, 
lim u(x, t) = 0, t > 0, 

x— >±oo 

where a is a nonzero constant, g and q are arbitrary smooth nonvanishing functions and n ^ 0, 1. 

We have derived the Lie symmetries for the variable coefficient equation ([TJ and now we 
examine which of these symmetries leave the initial and boundary conditions of the problem (I16p 
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invariant. The procedure starts by assuming a general symmetry of the form 



(17) 



where n is the number of Lie point symmetries of the given partial differential equation and the 
ai, i = 1, . . . , n, are constants to be determined. 

Lie symmetries for equation (TjQ) appear in Table 1. In case 2, for which g(t) = et p , the 
generator (I17p takes the form 



T = a\d x + a 2 



2td t + (p + l)xd x + 



1 



ud u 



n — 1 

Application of V to the first boundary condition which is written as 
x = and u(t, 0) = q(t) 

gives 



a\ = and a 2 — 2£— H g = 0. 

at n — 1 

For nonzero «2j we have 

where 7 is an arbitrary nonzero constant. It can be shown that symmetry, T, with a\ = 
leaves invariant the second boundary condition and the initial condition. Hence the admitted 
Lie symmetry can be used to reduce the initial- and boundary- value problem (|16p to a problem 
with the governing equation being an ordinary differential equations. In fact the Lie symmetry 
ltdt + (p + l)xd x + lf^jud u produces the transformation 

u = t 2n ~ 2 (p{r]), 7] = xt 2 



which reduces (1161) into 



(18) 

(19) 
(20) 



2e0" + (p + 1)#' - 2a(0 n )' - ^— U = 0, 

n — 1 

0(0) = 7, 
lim 0(77) = 0. 

77— >±oo 

Let p = (2 — n)/n. Then (|18p takes the form 

^" + I(#' + 0)-a(0")' = O 

and can be integrated once. Then we get ecf>' + \ f]4> — o,cf) n + c = 0. When we set c = 0, we are 
able to find a particular solution of the form 



e 2en 



^'C + il^^erff^ 



1 

l-n 



4e 



(21) 



where a = 2ye(n— TJ/n. This is the real- valued solution of the BVP (|18p - ([20p when C = 7 1 
and the condition en > holds. 
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Therefore we obtain the solution of the BVP 



ut + a(u n ) x = et n u xx , — oo < x < +00, t > 0, 
u(0, x) = 0, — 00 < x < +00, 
u(i,0)=7i~n, t>0 
lim u(x, t) = 0, t > 0, 

x— >±oo 



where e > 0, n > 1, in the form 



i « exp 



1 9 2 
2en 



n a 



27r erf 



a 



V2a 1 

XI n 

4e 



1 

1— n 



(22) 



where a = 2^f e(n — l)/n. This solution was found in pQ, but in our opinion, by the use of 
a more complicated method. 

The above procedure can be applied to the remaining cases that appear in Table 2. If we 
omit case 9 which is the well-known Burgers Equation, only in case 7 there exist a Lie symmetry 
that leaves the boundary and initial conditions invariant. However, the results for this case can 
be obtained from the above by setting n = 2. 
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